OG 3.3.25 



Cosmic Ray Acceleration by Magnetic Traps 

V.N.Zirakashvili 1 

1 Institute for Terrestrial Magnetism, Ionosphere and Radiowave Propagation, Russian Academy of 
Sciences (IZMIRAN), 142190, Troitsk, Moscow Region, Russia 

Abstract 

Cosmic ray acceleration in turbulent interstellar medium is considered. Turbulence is treated as 
ensemble of moving magnetic traps. We derive equations for particle momentum distribution function 
that describes acceleration of particles in this case. Rate of acceleration calculated is estimated for our 
' Galaxy and compared with ones given by other acceleration mechanisms. 

° : 

^ ■ 1 Introduction. 

Several types of distributed cosmic ray acceleration were investigated by many authors: resonant 
stochastic acceleration (Tverskoi,1968), stochastic acceleration by large scale turbulence (Bykov & 
Toptygin, 1979), friction acceleration (Berezhko, 1981, Earl, Jokipii,& Morfill, 1988, Webb, 1989) and 
in earlier papers of Fermi (Fermi, 1949) and Alfven (Alfven & Falthammar, 1963). The last author 
, introduced so called "magnetic pumping". For such an acceleration mechanism particle changes its 
perpendicular momentum in varying in time magnetic field, the gain of energy being redistributed due 
to scattering. In this paper we shall introduce further development of "magnetic pumping" and apply 
, it for cosmic ray acceleration by magnetohydrodynamic (MHD) turbulence. 

In our Galaxy interstellar MHD turbulence is mainly created by supernova (SN) explosions. There 
exists turbulent transfer of kinetic and magnetic energies from the main scale L ~ 100 pc to smaller 
scales (Ruzmaikin, Sokolov,& Shukurov, 1988). In this paper we shall treat this turbulence as large 
\ amount of moving magnetic traps. Energetic particles are trapped inside each trap and change energy 
Qh' due to trap deformation. Weak scattering slowly changes pitch-angle of the particle, determines escape 
q ■ time and hence rate of acceleration. 

h ; 

2 Acceleration of particles in magnetic traps. 

Assuming that particles gyroradii are much smaller than magnetic trap size we will use drift equation 
. £h ' with scattering (Chandrasekhar, 1960) 
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Here / is cosmic ray momentum distribution function that is averaged on gyroperiod, p is particle 
momentum, b = B/i? is the unit vector along the magnetic field B, fj, = pB jpB is the cosine of 
pitch angle of the particle, = c[E x B]/S 2 is the electric drift velocity. For frozen magnetic field 
it is simply perpendicular to B component of medium velocity u. The right hand side of equation (1) 
describes scattering of particles by small scale magnetic inhomogeneities moving with medium velocity 
u. Isotropic function f a that is independent on time and space coordinates is exact solution of Eq. 
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(1) for u = 0. Small deviation Sf of cosmic ray momentum distribution function arises for nonzero u. 
Linearizing Eq.(l) we have 
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Left hand side of Eq.(2) describes scattering and movement of particle, the third term being correspond 
to magnetic mirroring. The terms proportional to uu are neglected since we are going to find solution 
of this equation in the limit of small scattering frequency when the particle free path length A = vjv 
is much greater than the trap size I. In this case it is convenient to use new variable q = (1 — fi 2 )/B 
instead of \i. Let Sf = Sf + for /i > and Sf = Sf~ for \i < 0, then assuming quasistationarity of 
Eq.(2): 
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Here s is the coordinate along the magnetic field line. Our aim is to perform averaging of Eq.(l) on 
ensemble of realizations of random magnetic field and medium velocity. Since mirroring should be 
included the problem is nonlinear. It will be simplified if one assumes that the medium velocity scale is 
larger than the magnetic field scale and will treat them as statistically independent quantities. We will 
also assume that medium is filled by many identical magnetic traps, and magnetic field of each trap 
varies from B m \ n up to B max . The same treatment was used for investigation of influence of magnetic 
traps on the second harmonics of cosmic ray anisotropy (Klepach & Ptuskin, 1995). 
For trapped particles l/B max < q < 1/-B m in- These particles perform oscillations inside a trap and 
slowly diffuse on q. Hence function 5/^ for these particles is approximately independent on s and can 
be found from the equation: 
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Here s\(q) and S2(q) are roots of equation B(s)q 
following equation: 
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1. Performing one integration on q one obtains the 
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Substituting the solution of this equation into Eq.(l), performing ensemble averaging and integration 
on \i one can find momentum diffusion coefficient (energy changes of trapped particles are taken into 
account only): 
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Here w is the volume factor of traps and (...) means ensemble averaging on u. Expression (6) is 
rather cumbersome. In order to simplify it we assume that the magnetic field disturbances are small: 
5B/B « 1. In this case trapped particles have pitch angles close to 7r/2, that is fi « 1 and one can 
let qB(s) = 1 in square brackets of last integral. We will also neglect by second term in this integral 
in comparison with first one. Assuming also that medium velocity scale is larger than trap size I and 
using //-independent scattering frequency one can obtain 

D„ = ((Vu ± ) 2 ) £ / ds X l ( 1 - ) (7) 

The deviation of momentum distribution function is given by the expression: 

max ~~ B{s) 2 \ a ( -^max ~ B{s) 




It "jo " "_l / "max -— v" ) 2\ a I "max "\"J 2 \ /o\ 

S!= ^— {—B^ "A" "I (8) 



-"max 



where 9(x) is step function and uj_ is perpendicular to the mean field component of the medium 
velocity. 

The expression (7) can be readily obtained from the following estimation. For trapped particles 
< [i*. If trap is compressed or expands in perpendicular to the mean field direction particle 
changes its momentum with the rate ^| ~ pVu±. The escape time is At ~ fj£/is. Taking into account 
that trapped particles are only wfj,* part of all particles one can find that 

If magnetic field disturbances have component parallel to the mean field (magnetosonic wave magnetic 
disturbances) then //* ~ ^JbBjB and we obtain expression (7). 

3 Discussion. 

It should be mentioned that expression (9) has the same form as the one for the friction acceleration 
(Berezhko, 1981, Earl, Jokipii & Morfill, 1988, Webb, 1989). The friction acceleration was derived for 
the case when medium velocity scale is larger than particle free path length. The acceleration considered 
here works for opposite relation I << A and thus can be stronger for sufficiently small I. The same is 
concerned for comparison with stochastic Fermi acceleration because it is of the same order as friction 
acceleration for I ~ A. The rate of resonance acceleration is roughly (cf. Berezinsky et al.,1990) 

^~4(§) 2 ^ 2 (io) 

where f2 is particle gyrofrequency, v a is Alfven velocity, and SB is the value of random magnetic field 
in the scale of the order of particle gyroradius r g = v/Cl. In order to compare estimations (9) and (10) 
one should take into account that u ~ v a 5B/B, //* ~ ^SB/B, w ~ 1 and take the smallest possible 
I ~ r g when the drift approximation is marginally valid. Taking for the scattering frequency value that 
is determined by resonance scattering u ~ £1 (5B/B) 2 (cf. Berezinsky et al., 1990) and assuming that 
velocity scale is of the order of the trap size one can find that 
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Therefore resonant acceleration is a factor (5B/B) ' less effective than considered here trap accelera- 
tion. If scattering frequency drops at small fi the difference will be larger in favor of trap acceleration. 



It is well known now that resonant acceleration is enough for weak reacceleration of galactic cosmic rays 
with energies about lGeV/nucleon in order to explain mean thickness dependence at these energies 
(Seo & Ptuskin, 1994). Problems can arise for more powerful acceleration. It seems that magnetosonic 
waves with lengths of the order of 1 GeV proton gyroradius cannot exist in the interstellar medium 
because strong Landau damping (acceleration rate estimation is obtained for such waves). 
In our Galaxy compressible turbulence exists in the main scale L ~ lOOpc and possibly at smaller 
scales due to nonlinear cascade transfer. The minimum scale of compressible turbulence is determined 
by dissipative processes and is estimated as L m \ n ~ 0.01 — lpc (Ruzmaikin, Sokolov, & Shukurov, 1988). 
It should note that this minimum scale can be also determined by dissipation related with the accel- 
eration mechanism considered. 

The acceleration rate (9) is determined by values of velocity and magnetic perturbations for this 
minimum scale. In this case mechanism of acceleration considered has very attractive feature: for 
ultrarelativistic particles the rate of acceleration is proportional to space diffusion coefficient. Such a 
feature gives a possibility of formation of power-low cosmic ray spectrum. This dependence is valid up 
to the energies when particle gyroradius is approximately equal to minimum scale of turbulence and is 
estimated as Z ■ lO^eV where Z is the particle charge, -/^min ~ O.lpc and B ~ 10~ 5 Gs was assumed. 
The ratio of diffusive exit time to acceleration time can be estimated as 



Here Hh is the galactic halo height, h acc is the height of acceleration region, SB m { n is the value of mag- 
netic field perturbation for minimum turbulence scale L m [ n . This ratio is about 0.1 for L m [ n ~ O.lpc, 
Hh ~ Wkpc, h acc ~ 3kpc, v a ~ lOOkm/s, 5B m \ n /B ~ 1/10. 

For rough comparison with acceleration on SN shocks one should note that both mechanisms use such 
a source of energy as mechanical energy of SN explosion. Since only about 5 percent of initial mechan- 
ical energy of SN explosion goes into mechanical energy of interstellar medium (Ruzmaikin, Sokolov 
& Shukurov, 1988) one can expect that if more than 5 percent of the initial mechanical energy of SNs 
goes into accelerated at SN's shock particles, then acceleration mechanism considered here is negligible 
in comparison with SN shock acceleration. In the opposite case it can be important for our Galaxy. 
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